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I. INTRODUCTION 

Two-component fermionic quantum gases belong to 
the most versatile systems that can be used to study 
the many-body physics of fermions [TJ [5] . The simplest 
Lagrangian to describe such a system models the s-wave 
interaction between the constituents by a contact inter- 
action (see, for example, Ref. [3]): 

C = V4 (idt + A^Vj^t- (!) 

where a summation over spin indices a =t, i is implied. 
In experiments with cold fermionic atoms, it is possible 
to confine the atomic cloud in two- and three-dimensional 
traps and also to tune the interaction strength to a regime 
where the components are strongly interacting. The the- 
oretical analysis of a Fermi gas in a strongly interacting 
limit is intricate since many theoretical methods rely on 
a perturbative expansion or neglect fluctuations and are 
thus not applicable. 

Over the past few years, a lot of work has been stimu- 
lated by the observation that some properties of a Fermi 
gas can be described by universal relations. They are 
called universal because they hold over a wide range of 
the system's parameters such as the temperature, the 
density, the precise form or the strength of the interac- 
tion, and so on. A first set of such relations was pre- 
sented by Tan [1HB], who among other things derived 
expressions for the momentum distribution, the energy, 
and various other thermodynamic quantities, and are 
hence often called Tan relations. The scale of the Tan 
relations is set by a quantity known as the contact C, 
which is determined by the many-particle properties of 
the system. It is defined as C = Jd d R(C(R)), where 

C(R) = m 2 A 2 V ; |''/']'0^'0^(R') is the contact density opera- 
tor and d = 2, 3 denotes the number of space dimensions. 
An operational definition of the contact can be given as 
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the magnitude of the large-momentum tail of the gas: 
C = lim/ i ;_ s . 00 n(k)k i . Apart from the dependence on the 
contact, the universal relations can be inferred from the 
few-body physics and can be calculated in closed analytic 
form in many cases. Only the value of the contact has 
to be taken from numerical calculations or experiments. 
Some of the Tan relations have been tested experimen- 
tally, and good agreement with the theoretical predic- 
tions has been found [7j. Since their proposal, many other 
exact relations for fermionic as well as bosonic systems 
have been derived (cf., for example, Refs. [5HT7]). An 
overview of various universal relations and recent exper- 
imental developments can be found in Ref. [18] . 

A field-theoretical derivation of universal relations was 
put forward by Braaten, Kang, and Platter [19, 20 , who 
used the operator-product expansion (OPE) to rederive 
the Tan relations [IHB]. The OPE expresses the product 
of two local operators and at different points 

(k) 

in space and time in terms of local operators O a , where 
a = (cti, ot2, ■ ■ ■) denotes a set of indices [2"T] : 

Of (R + r , T + t)Of (R, T) = Y, c ?i T ) ■ 

k,a 

(2) 

Note that, on the right-hand side, the dependence on the 
separation (r, t) is carried by the Wilson coefficients c^ a , 
which can be calculated exactly for simple operators. If 
an observable can be expressed as the expectation value 
of an operator product involving operators at different 
space-time points, we can use the OPE to determine its 
leading-order behavior for small r and t (or for large val- 
ues of momentum q and energy uj for the Fourier trans- 
form). The advantage of the OPE is that it provides a 
powerful framework to derive additional universal rela- 
tions. This was exploited by the authors of [2U] to derive 
universal results for radio-frequency spectroscopy mea- 
surements |llj . Son and Thompson studied the asymp- 
totic behavior of the dynamic structure factor S*(q, uj) at 
unitarity at large q and uj as well as the large-frequency 
tail of S [12]. The result for the large- frequency tail in 
[T2"] was augmented by Goldberger and Rothstein, who 
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calculated the leading-order correction due to finite scat- 
tering length (55]. Some of the above results have also 
been obtained by other authors using different methods 
[TH - fTt)] . While most of the above work was carried out 
for a three-dimensional system at or close to unitarity, it 
is possible to derive universal relations for any strength 
of the interaction or any other dimension in the same 
way. 

In this work, we apply the OPE to derive the current 
response function of a Fermi gas in two and in three di- 
mensions. It is defined as 



Gf/ (x-y,t) =ie(t)<[tf(x,t),# (y,0)]), (3) 



where jf(x,t) = -(i/2m)(ipldiip a ) = -(i/2m)(ipldiip a 
—di^l-ipa) is the current operator. The Fourier transform 
of Eq. @ can be split into a transverse and a longitudinal 
part: 



G^(q, W ) = Gr (q,w) 



9iQj 



G° T ° (q,w) 



The longitudinal part is related to the superfluid fraction 
and the transverse part to the normal fraction of the gas. 
Thus, the current response function is a central quantity 
in the study of superfluids. The longitudinal component 
is related to the structure factor and can be determined 
with methods that measure the density response, such as 
two-photon Bragg spectroscopy [25] . 

Recently, there have been several proposals to probe 
the transverse component of the current response using 
artificial gauge fields [2"1H2"B] . Artificial gauge fields are 
induced by either rotating the atomic cloud or placing 
it in a specific laser setup. In the latter case, incident 
laser beams couple the ground state of the atoms to an 
excited state. The effective Lagrangian of such a system 
describes charged particles coupled to a U(l) gauge field. 
The amplitude of this gauge field typically depends on the 
spatial or temporal variation of the Rabi frequency or the 
detuning of the laser from the resonance frequency. For 
example, a local detuning of the laser field could induce 
a rapidly oscillating gauge field. For further references, 
see [571 HH] • One way to measure the imaginary part of 
the current response is to determine how much energy W 
is absorbed by the system if the gauge field is switched 
on for a short time: 



W 



duj 
, 2^ 

x 4i/(q,w 



d d q 

(2~7rp 



Al(-q, -u)ulmGff (q,w) 



(5) 



For a two-dimensional gas, there also exists a proposal 
to determine the current response optically, by measur- 
ing the phase shift of a laser beam that crosses the atomic 
cloud [2S] . While artificial gauge fields have only been re- 
alized for Bose gases so far there is no reason why this 
should not also be possible for fermionic systems in the 
near future. It is in anticipation of such experiments 



that we calculate the current response. In particular, 
since it should be possible to engineer gauge fields that 
couple differently to each spin species, we consider both 
the spin-symmetric and antisymmetric current response 
G'j" — Gfj 7 ± G?j~" , where a summation over spin in- 
dices is implied [corresponding to gauge fields that cou- 
ple symmetrically and antisymmetrically to the current: 
A-(j t ±jJ. 

This paper is structured as follows: In Sec. |TTJ we de- 
scribe how to calculate the Wilson coefficients of the cur- 
rent correlator. This is done using a matching proce- 
dure where the matrix element of Eq. ([2]) is evaluated 
with respect to selected few-body states. We choose a 
one-particle state to determine the Wilson coefficients of 
the bilinear operators (Sec. II B[ ) and then a two-particle 
state for the contact operator (Sec. II C). Section III 
presents the results of this calculation. Expressions for 
both the spin-symmetric and the antisymmetric current 
response are presented and discussed in Sec. |III A| The 
U(l) symmetry of the system implies a Ward identity 
that relates the longitudinal part of the current correlator 
to the dynamic structure factor. We use this to calculate 
the asymptotic form of the dynamic structure factor in 
two and in three dimensions for any value of the scatter- 
ing length in Sec. IIIB In [11], Taylor and Randeria de- 
rived expressions for the large-frequency tail of the spec- 
tral shear and bulk viscosity rj(ui) and This paper 
reproduces their results using a field-theoretical frame- 
work and generalizes them to finite scattering length. 
We are also able to provide the viscosity tails for a two- 
dimensional Fermi gas. We conclude with a summary in 
Sec. |IV| The paper contains an appendix that presents 
details of the calculations in Sec. HU 



II. OPERATOR PRODUCT EXPANSION 



While we would like to obtain the OPE for the re- 
tarded response function G°° (q, w), it is more conve- 
nient to calculate the time-ordered response, which co- 
incides with the retarded response up to a correction of 
order 0(2e~^) [25]: 



G°f{q,uj) = i J dte i < w+ie >* J d d r J d d R e - <q " 



x (Tjf (R+ ~,t) jf (R - |, 0)} + 0(2e~^). 



(6) 



T denotes the time-ordered product. The it term in the 
Fourier transform has been included to shift the poles in 
the resulting expressions of the response functions [25] . 
Using the time-ordered response allows us to apply stan- 
dard diagrammatic methods to calculate the Wilson co- 
efficients. In this section, we derive the OPE of the op- 



3 



erator 



Acra 



d d r I d d Re lwt ~ iqr 



(q,w) = J dt 
xTjf (R + |,t) #' (R-|,0 



(7) 



The OPE is an operator identity, that is, it holds if 
we take the expectation value of Eq. ^ with respect 
to any state. Hence, in order to determine the Wilson 
coefficient of an operator 0^ k \ we can choose a sim- 
ple few-body state for which (O^) =^ and match the 
expectation values on both sides of Eq. 0. Inserting 
Eq. in Eq. yields 



AtJ<J 



E 



i 



„(*0 



w A fc /2-d/2 



2mcj ' 



y d d R0^(R), 



(8) 



where denotes the scaling dimension of the operator 
o£ fc) as defined by oi fc) (Ax,A 2 £) = X-^O^ (x,t). We 
suppress the dependence of the Wilson coefficients on the 
spin indices on the right-hand side of Eq. 0. The op- 
erators Oq' can be composed of the field operators ip a 
and their adjoints as well as derivative operators act- 
ing on those fields. We will only consider operators that 
conserve the particle number, i.e., contain the same num- 
ber of fields and their adjoints, and that are symmetric 
in the spin indices. Equation shows that the high- 
frequency behavior of the current response is governed 
by the operators with lowest scaling dimension. Those 
operators are 



n = V4Vv (A„ = 


d), 


(9) 


J = -^4VVV (Aj = 


d+1), 


(10) 


Off = ^£ (A <3, 


= d + 2), 


(11) 




= d + 2), 


(12) 


and the contact density operator 






n ty 4. A 

C = m X ip^xp^ip^ip^ with 


A c =4. 


(13) 


<-> 

We use the shorthand ifp d\f = ij^dip 







We proceed as follows: First, we summarize the Feyn- 
man rules of the theory in Eq. and match the bi- 



linear operators 0-( 12 1 using a one-particle state (for 
which the matrix element of the contact density van- 
ishes). Then, we determine the Wilson coefficient of 
the contact using a two-particle state and the results of 
Sec. ITTBl 



A. Feynman rules 

The matrix element of an operator can be expressed 
pictographically in terms of Feynman diagrams. This 
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FIG. 1: (a) Feynman diagram for the two-particle scattering 
amplitude, (b) Current insertion on a fermion line. 



section lists the Feynman rules that are necessary to eval- 
uate those diagrams. 

Every line is assigned an energy and a momentum 
where we impose momentum and energy conservation 
at each vertex. The external energy and momen- 
tum are constrained by the state. After applying the 
Feynman rules, we integrate over each undetermined 
loop momentum k and energy fc using the measure 
J d d 1idko/(2ir) d+1 . Each internal line of a particle with 
energy ko and momentum k contributes a factor i / (kp — 
Ek + ie) to the integrand. The dashed blob [Fig. 1(a)] 
denotes the two-particle scattering amplitude and gives 
a term i^4(fco,k), where ko and k are the sum of the 
energies and the momenta of the incoming fermions, re- 
spectively. We have 



A(k, k ) 



4?r 



1 



2 ~ a 3D + y/~ m (ko — k 2 /4m) — ie 
in three dimensions (3D) and 

An i ^ 2n 1 
AOs., k ) 



m lna 2 -Dy / -m(ko - k 2 /4m) 



■ ie 



(14) 



(15) 



in two dimensions (2D). In the following, wc abbreviate 
A(0, 0) = A. The strength of the interaction is given by 
the 2D and 3D scattering lengths a^o and a^D, respec- 
tively. ,4(fco,k) can be obtained by summing all ladder 
diagrams that contribute to the two-atom Green's func- 
tion. The bare coupling constant A is related to A by the 
rcnormalization condition 



A~ 



1 

A 



d d k 1 
(2^2^- 



(16) 



This condition ensures that the two-particle scattering 
amplitude reproduces the leading order in the effective 
range expansion of the quantum-mechanical s-wave scat- 
tering phase shift, fccot<5(fc) = — in 3D [3] [30] and 
cotd(k) = (l/7r)lnA; 2 a 2 D in 2D J32GD2], respectively. 

We denote the insertion of a current operator by a 
crossed circle attached to a wiggly line. A current inser- 
tion with energy u) and momentum q that is attached to a 
fermion line with energy k and momentum k is depicted 
in Fig. 1(b) It contributes a factor 



1 

2m 



(2fci + qi) 



(17) 



D. On 



D. On 



FIG. 2: Diagrams that contribute to the expectation value 
of Alf (q, oj) [Eq. ^] between a one-particle state with mo- 
mentum p and energy po- 



FIG. 3: Insertion of a bilinear operator on an internal line. 
The insertion on an external line matches the insertion of two 
current operators on the same line and need not be considered. 



to the integrand. The Feynman rules of the opera- 
tors ([9j)-([T2]) are listed in Eqs. ([l9|-([22]). 



B. One-particle state 

In order to determine the Wilson coefficients of the bi- 
linear operators, we evaluate the matrix element of Af? 
[Eq. |7])] with respect to a one-particle state with momen- 
tum p and energy po. Note that we do not enforce the 
on-shell condition po = p 2 /2m. The relevant Feynman 
diagrams are depicted in Fig. [2j Their value is: 



(p,Pol^l 7 T (q,^)|p,Po) = 



i/Am 2 

J po+oj - E. 

x (2 Pl + qi)(2pj + qj) + (w, q 



p+q 

w,-q), 



(18) 



where E p = p 2 /2m. The same holds for Ajj. The ex- 
pectation values of the bilinear operators are 

(p,Po\n\p,Pa) = 1, (19) 

(p,Po\jcx\p,Po) = Pa/m, (20) 

(p,Po|oi 3 j|p,?3o) = -2p aPp /m and (21) 

(p,p o |0 (4) |p,Po) =Po-E p . (22) 



They match the first terms in the expansion of Eq. ( 18 ) 

5 yields t 



in terms of po and p. This yields the Wilson coefficients 
i/Am 



c J - = 



moj — q 2 /2 
i/4 



(23) 



(m W -q 2 /2) 2 " 9afi * ' mw-q 2 /2 
x (S ia qj + 5j a qi) - (w -> -w), 
2 - 1 12 



i/2 



,(3) 



I X 

a f) d *f3 



(24) 



-ijaP 



-iq 



{mu-q 2 /2y {muj~q 2 /2) 2 
i/2 



and 



J 4 ) 



mw - q 2 /2 *" 
-i/4 



% + (uj -t -uS) 



{muj — q 2 /2)' 



qtq-i + (u) -> -w). 



(25) 



(26) 



For an isotropic system, we have \0 K a L) = (2/md) 
(■0j.V 2, ( j)^Q0, hence, only the contraction c^^<J Q| s will 



be relevant in the following. Furthermore, we have 
(j) = and can apply the Heisenberg equation of motion 
of the operator O^: (OW) = (2/to 2 A)(C). The OPE 
allows us to express the expectation value of Eq. ^ as 
follows: 



(Aff (q, w)> = c" (q, w) N + & (q, w) H 



+ cf ? (q,u,)C + 



1 



(27) 



Here, N = J d rf R(n(R)) and H = j d d H (%(R)), where 
% denotes the Hamiltonian density 

W(R) - ^vv4vv CT (R) + AVM^MW, ( 28 ) 



which has the Wilson coefficient 



<^(q,w) = -j8 a pc$ afj (ti,w). 



(29) 



From Eq. (27 1, we see that the terms proportional to 
N and H only contribute a term ~ 5(muj — q 2 /2) to 
the imaginary part of Gf? . Away from this one-particle 
peak, the leading order is proportional to the contact. 
In the next section, we proceed to calculate the Wilson 
coefficient cf (q, w) of the contact operator. 



C. Two-particle state 

To determine the Wilson coefficient of the contact op- 
erator, we choose a two-particle state where both par- 
ticles have zero energy and momentum. This way, the 
matrix element of other two-body operators with higher 
scaling dimension than the contact, which involve addi- 
tional derivatives of the field, are zero. We only have 
to subtract the matrix elements of the one-particle op- 
erators with respect to the two-particle state. We do, 
however, introduce additional infrared divergences when 
calculating single diagrams, which we need to subtract 
before we can evaluate the finite part. In particular, in 
two dimensions, this makes the calculations somewhat 
lengthy. Of course, the overall result is finite and the 
divergent parts cancel when summing all diagrams. 

We start by calculating the insertion of the bilinear 
operators on an internal line as shown in Fig. [3] To 
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obtain a finite expression, we introduce an IR and a UV 
cutoff 7/ and A, respectively: 



(0,t;0,||n|0,t;0,|)| Fig[ 3]= (n) 

i/4irr] d = 3 
— 1/2tt77 2 d= 2 ' 



mM x 



(30) 



(0,t;(U|H|0,t;(U>U 



m 2 A 2 



l 

2tt 



3D 



/4tt 



(31) 



and 



(0,t;0,||O (4) |0,t;0,|) 

2m 2 A 2 



= (0< 4 >) 

Fig-! 

A/2vr 2 = 3 

a2 1 rf = 2 



The matrix element of the contact operator is 

(0, t; o, 4 I C | o, t; o, 4) = (C) = m 2 ^l 2 




(32) 



(33) 



both in two and in three space dimensions. To keep our 
notation concise, we abbreviate the matrix elements of 
the various operators as defined in the previous equa- 
tions. 

The relevant diagrams that contribute to the matrix 
elements of Eq. ([7]) are shown in Fig. [4 The insertion 
of two current operators on an external line has already 
been considered in the previous section and is matched by 
bilinear insertion on external lines. Figures [4(i)1and[4(b)| 
respectively, contribute to only the diagonal and off- 
diagonal components of (A^f ) in spin space, whereas 
Fig. |4(c) contributes to all components. The matrix ele- 
ments of the bilinear operators are diagonal in spin space. 
Thus, they cancel the divergent parts in Fig. 4(a)| Fig- 
ures |4(b)| and |4(c) are separately fini te. In the f ollo wing, 
we denote the contribution of Figs. 4(a) |4(b) and |4(c) 



to the Wilson coefficient of the contact by cjj, C-^ and 

C\j , respectively, where we absorb the contribution of 

the bilinear matrix elements into the definition of C. 
In short, we have 



(i) 



4 = C% } + C\f ^ + Cf ■ (34) 

Using the results of this section, we can give an explicit 
expression for the coefficient : 



« 2m 4 J (2n) d 4E 2 u - £ k - £ k4 



is 



X (2k, + qi )(2k 3 + q 3 ) - c?M - c « (q , 



.(4)/ 



; 2 A 



(C) 



+ (35) 




w, q 




(c) 

FIG. 4: Diagrams that contribute to the matrix elements of 



A™ (q, u) [Eq. |7l] between the two-particle state. 



The bilinear contributions cancel the IR divergences of 
Fig." 



4(a) and render the expression (35 1 finite. In 3D 



and 2D, Fig. 4(a) has a power-law infrared divergence 



that is matched by the matrix element of the number op- 
erator n. In 2D, there are additional logarithmic infrared 
divergences in both the transverse as well as the longitu- 
dinal components. They are canceled by the expectation 

values of the operators % and . 

The sum of the diagrams in Fig. |4(b)| (divided by the 
matrix element of C) is equal to 



C 



(2) 



d d k 1 



1 



x (2k t 
+ (u>- 



(2ir) d E^lo- 
qi)(2kj + qj) + 



2iA~ x 



■ IE 



(36) 



In two dimensions, the first term contains a logarithmic 
infrared-divergent part. It is canceled by the infrared 
divergence in the two-particle scattering amplitude A [cf. 
Eq. pj]. 

The last contribution to c£ comes from the three dia- 



grams in Fig. 4(c) They only contribute to the longitu- 
dinal part of the response (cf . the Appendix) : 



C 



(3) 



L4(w,q) 



u) - E a 



A~ 



d d k 

(27r) d 



G 



1 



1 



+ (cj -> -w). 



■ 2£ 



(2kj + ft) 



x [t j] 
(37) 



As before, a logarithmic infrared divergence in 2D is can- 
celed between the two terms and the expression is mani- 
festly finite. 

It is possible but cumbersome to evaluate the expres- 
sions (35), (36), and (37) in closed analytical form in 



terms of elementary functions. For details of this cal- 
culation we refer to the Appendix. The results of the 
computation are presented in the next section. 



III. RESULTS 

The first part of this section lists and discusses the 
current response functions of a Fermi gas in two and 
in three dimensions for large values of the arguments 
q and cj, where x = q 2 /2muj and a~ x j \Jmu are kept 
fixed. We consider both the response to a gauge field 
that couples symmetrically as well as antisymmetrically 
to the different fermion species. The latter corresponds to 
probing the response of the spin-antisymmetric current 
j a = y — y. We denote the spin-symmetric and anti- 
symmetric response function by a superscript s and a, 
respectively. In the remainder of this section, we use the 
results of Sec. |III A| to calculate the dynamic structure 
factor in 3D and 2D. We also determine the asymptotic 
form of the spectral viscosities introduced in [T3] for large 
values of the frequency, and comment on the derivation 
of sum rules using the OPE results. 



A. Response functions 



As discussed in Sec. II C cj^ [Eq. (35)] contributes to 
the diagonal and - [Eq. (36)] to the off-diagonal com- 



ponents of the con tac t's Wilson coefficient in spin space, 



whereas C z - 3 ^ [Eq. (37)] is independent of the spin indices. 



Hence, we obtain the spin-symmetric and antisymmetric 
response by summing those contributions as follows: 



lmG s T/L {q,uj) 



J T/L 

~c 



= P i / Im 



id) 1 + iC)V + iC. 



(2) 



,(?•) 



and 



ImG^ /L (q,w) 
C 



i^ /T Im 



iC 



(i) 



iC. 



(2) 



(38) 



(39) 



where Pfi = q t qjjq 2 and = (% - q^ / q 2 ) / (d - 1) 
are the projectors onto the longitudinal and the trans- 
verse parts, respectively. Remember that we have to 
shift cj — > cj + is when evaluating the Wilson coefficients 
[cf. the definitions in Eqs. j6l) and Q]. As discussed 

(3) 

in the Appendix, of the above terms, only ImiCj-- de- 
pends on the scattering length away from the one-particle 
peak. Barring the dependence of the contact on a -1 , this 
means that, except for the spin-symmetric longitudinal 
response function, all response functions are independent 
of the scattering length in the asymptotic regime. 

The imaginary part of the response functions presented 
in the next two sections vanishes for x > 2. As al- 
ready pointed out in [T^] , this represents the two-particle 
threshold. This lower bound lj > q 2 /Am follows from 
energy and momentum conservation if we consider the 
absorption of a probe with (large) momentum q and en- 
ergy u) by two particles in the medium. 



3D 



The contributions of Eqs. (35)-(37) to Eqs. (38) 



and ( 39 ) can be inferred from the results in the Ap- 
pendix. We start by listing the response functions for 
the three-dimensional case. The spin-symmetric and an- 
tisymmetric current responses are 

Im Gy(q, cj) 



C Aitrnq 2 
2{muj) 2 + (mu - q 2 /2) 2 
4:Trm 2 ujq 3 



\J muj — q 2 /A 

muj + q\J moj — q 2 /4 



In- 



\mcj - 



q 2 /2\ 



(40) 



and 
ImG^(q,w) 

2(muj) 2 



Airmq' 



\/muj — q 2 /A 



[mcj — q 2 /2) 2 ^ muj + q^Jmcj — q 2 /A 



Anm 2 ujq 



\muj - q 2 /2\ 



(41) 



Figure [5] shows plots of the transverse response func- 
tion as a function of th e scal ing variable x = q 2 /2muj 
[Fig. |5(a)| and a; -1 [Fig. 5(b) . They can be interpreted 



as plots at fixed energy and momentum, respectively. As 
mentioned earlier, the transverse components do not de- 
pend on the scattering length. 

Similarly, we derive the longitudinal response functions 



lmG a L (q,cj) 
C 



q 2 /A 



2irq 2 (mcj — q 2 /2) 2 2Trq : 



In 



mcj - 



mu — q 2 1 A 
V2j 



(42) 
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m Im[G7-(-v)] 





FIG. 5: 3D: Transverse response function at (a) fixed energy and (b) fixed momentum as a function of the variable x = q 2 /2rnui 
and respectively. The continuous blue line denotes the spin-symmetric and the red dashed line the spin-antisymmetric 
response. 



2 lm[G](x)] r r „ . _k, 

mqlm[Gi(x )] 




(a) (b) 

FIG. 6: 3D: Longitudinal response function at (a) fixed energy and (b) fixed momentum as a function of the variable x = q 2 /2rnui 
and x~ , respectively. The continuous blue line denotes the spin-symmetric and the red dashed line the spin-antisymmetric 
response. The spin-symmetric response is evaluated at unitarity. 



and 



Im (q^q^) 
C 



2irq 2 (mu~q 2 /2) 2 



2a 



-l 

3D 



mu> — q 2 /2 



a 3D<? 



u) mui + fli/mw — o 2 /4 
In 



2irq 3 



\mu> — q 2 /2\ 



q a 3I 5 + raw - q 2 /4 \mu> - q 2 /2 
1 - 



2tt9 (q 2 /2 - muS] 



imua^/q ' 1 



"3D 



+ raw - <j 2 /4 \mui - q 2 /2 



a 3D q 



n<d{q 2 /2- 



muj 2 \J mui — q 2 /<i 
nq 4 + mui - q 2 /4 

At unitarity, this result simplifies considerably: 

ImG[(q,u) mui 2 \fm, 



In' 



muj + q\J mui — <? 2 /4 



\mui 



q 2 /2\ 



-Ti 2 Q{q 2 /2-mu) S j. 



(43) 



q 2 /4 



C 



2irq 2 (mw-q 2 /2) 2 
mui 2 1 



u mui + q\l mui — q 2 /4 
In 



q 2 /4 



2nq 3 \mu~q 2 /2\ 
2 mjjj + q\J mui — (? 2 /4 



In 



q 2 /2\ 



it 2 8 (q 2 j2 - mu)j . (44) 
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1.80 



(a) 



(b) 



FIG. 7: 3D: (a) Longitudinal response function as a function of x and inverse scattering length a 3 ^/y/muj. (b) Detail of (a) 
at unitarity. The peak at unitarity is due to the strong interaction. The thick blue line denotes the response at unitarity and 



corresponds to the thick blue line in Fig. 6(a) 



nrto Im[G7-(jc)] 
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FIG. 8: 2D: Transverse response function at (a) fixed energy and (b) fixed momentum as a function of the variable x = q 2 /2muj 
and x~ , respectively. The continuous blue line denotes the spin-symmetric and the red dashed line the spin-antisymmetric 
response. 



Figure [6] shows plots of t he lon gitudinal resp onse f unction 



2. 2D 



as a function of x [Fig. |6(a) and x 1 [Fig. 6(b) , where 



Im G* L (q, u) is evaluated at unitarity. In Fig.|7J we show For the tran sverse response function, we find the result 



plots of the longitudinal response function as a function 
of both x and the scattering length aJ-^/y/muj. As it has 
already been pointed out in |12j . the peak at x — 2 in 
Fig. 6(a) is due to the strong interactions between the 
spin sp ecies at unitarity. This is also nicely illustrated 
in Fig. 7(a) where we find a peak at = and a 
finite response otherwise. A detail of the peak is shown 



in Fig. [7(b)] 



Im Gy(q, (jS) 



8m 2 u) \muj — q 2 /2\ mq- 



8Mr 



for the spin-symmetric and 

ImG%(q,a;) _ 
C 



muj — q 2 /4 
2m 2 Lj\m,LO — q 2 /2\ 



(45) 



(46) 



for the spin- antisymmetric current. Bot h func tions are 
plott ed in Fig. [8j as a function of x [Fig. 8(a) and x" 1 
[Fig. [8(b)1 , as before. 
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1.5 



2.0 



FIG. 9: 2D: Longitudinal response function at (a) fixed energy and (b) fixed momentum as a function of the variable x = q 2 /2mw 
and x~ x , respectively. The continuous blue line denotes the spin-symmetric and the red dashed line the spin-antisymmetric 
response. The spin-symmetric response is evaluated at a2o^rnw — 5 and a2oq = 5. 




nrio lm[G L (x,a 2 1 D )] 

C 




FIG. 10: 2D: Longitudinal response function as a function of x and the 2D-scattering length a^/y/muj. There is a peak at 
( a 2D = 0,x = 2) due to the enhanced intera ction in this limit. The thick blue line denotes the response at a^/y/mto — 0.2 
and corresponds to the thick blue line in Fig. 9(a) 



As for the three-dimensional case, the spin-symmetric 
longitudinal response function depends on the two- 
dimensional scattering length 02D: 

ImGJ(q,w,a^) 



C 



ujq 



2mu> 



e 



<t_ 

\muj-q 2 /2\ 3 g 2 (mw-q 2 /2) 2 ° I 2 
1 4mw 2 



7T 2 + In 2 a 2 D (mu; - g 2 /4) g>w - g 2 /2) 2 



Amu 



In 



mw — q 2 /4 
| ma; — q 2 /2| 
1 2 

x lna2rj(mu; — q 2 /4) 



e 



(47) 



The spin-antisymmetric response takes a much simpler 
form 



ImG2(q,w) 
C 



9 2 /4 



2 |mw-g 2 /2| 3 ' 



(48) 



Figure [9] shows plots of t he lon gitudinal respo nse fu nction 



as a function of x [Fig. 9(a) and x [Fig. 9(b) . Fig- 



ure [10] shows a plot of the longitudinal response function 
as a function of x and the scattering length a^j y/muj. 

It is interesting to note that, for infinite scattering 
length = 0, there is also a peak at x = 2, de- 
spite the fact that the two-dimensional Fermi gas is not 
strongly interacting. This reflects the fact that the scat- 
tering between the particles is enhanced in this limit, as 
it can be seen from the two-particle scattering ampli- 
tude in Eq. ( 15 ). The large value of the scattering length 



10 



serves to balance the small energy above the two-particle 
threshold to — q 2 /4m. 

When comparing Figs. [4] and [6] to g] and [9j we see 
that the width of the single-particle peak is much larger 
in 2D than it is in 3D. In 3D, the decay of the Wilson 
coefficients away from the pole at x = 1 has a logarithmic 
and a power-law behavior ~ |1 — x\~ 2 for the transverse 
and the longitudinal responses, respectively, as opposed 
to higher power-law decays ~ |1 — and ~ |1 — a;|~ 3 
in 2D. This is to be expected since we integrate over one 
space component less in 2D. 



agrees. The authors of |Hj also calculated the correction 
to the tail at leading order in a^. Again, our result 
is in agreement with theirs. Our calculation generalizes 
the cited results to finite scattering length. For the two- 
dimensional Fermi gas, we find no dependence on the 
scattering length in the high-frequency tail 

lim lim 5 (2D) (q,w) = —± — r G. (53) 
u->oo q->o 87rmo; 4 

C. Spectral viscosities 



B. Dynamic structure factor 

The Ward identity that follows from the U(l) symme- 
try of the Lagrangian in Eq. ([I]) gives the first hydrody- 
namic equation describing the conservation of mass, 



5 t (n)+V-(j> = 0, 



(49) 



where j = jT + y. If we substitute the Fourier trans- 



form of Eq. ( 49 1 in the spectral decomposition for the 



imaginary part of the current response function 



ImGj (q,w) = 7r (l - e 



-0U,\ 



a,b 



-PE a 



;a|j i (q)|6)<6|^(-q)|a>%-S 6a ), 



(50) 



where E\, a = Ef, — E a and E a are the eigenvalues of the 
Hamiltonian of the system, we can relate the longitudi- 
nal part of the spin-symmetric response to the density 
response G p (where p = mn is the mass density) and 
deduce an expression for the dynamic structure factor 

1 Til 2 O 2 

S(q,w) = -ImG p (q,w) = — ^-Im GUq, u)) . (51) 

7T TTU!' ! 



Comparing this with Eqs. (43) and (47 1, we obtain the 



asymptotic behavior of the dynamic structure factor. In 
the three-dimensional case, the result coincides with the 
expression at unitarity 12] [cf. Eq. (|44[)]. The general 
behavior of the structure factor is the same as for the 
longitudinal response function, Figs. [7] and [TUJ and we 
refrain from plotting it again. 



When expanding the expression in Eqs. (43 1 and (47) 



to leading order in q, we get the large frequency tail of 
the dynamic structure factor for mq 2 <C to: 



lim lim S^f^id) 



1 



it 2 m 1 / 2 to 7 / 2 



4 
45 



30 



( a 3DA 



1 



G. (52) 



At unitarity, this agrees with the result by jl4j up to 
a factor of 2/3. The same frequency tail was derived 
by the authors of [12 and (22) . with whom our result 



The authors of [14 derived Kubo formulas for 
frequency-dependent generalizations of the bulk and the 
shear viscosity in three dimensions. The starting point 
of their derivation is the Euler equation, 



(54) 



where Ily is the stress tensor for a viscous fluid in d space 
dimensions: 



IF; 



pSi 



(55) 



■mj — yvij t pViVj „ v , 
p denotes the pressure, v the velocity, p the density, and 



diVj + djVi - ' v ) 



+ CMV-v) (56) 



is the viscous stress tensor [33] . Retracing the derivation 
in [M], it is not difficult to generalize their result to two 
space dimensions. We obtain the relations 



m 2 uj . 



Re r)(ui) = lim — ^— Im Gfi(q, uS) 



and 



g^O q 



, . Re r)(uj) +ReC(w) = lim — — 
a J <?->o q z 



(57) 



ImGi(q,w) 



(58) 

The leading order of r)(to) and C(w) at large frequency 
is proportional to the contact. Comparing this with 
Eqs. p0|), (|43|, ([45]), and ((47}, we obtain the large- 



frequency tails of the bulk and the shear viscosity: 

C 



and 



ReC(w) 



C 



in three dimensions, and 

Re 77(0;) 
and 

ReCH = 



C 
8mio 



(59) 

(60) 

(61) 
(62) 
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in two dimensions. The high-frequency tails have been 
considered before for a three-dimensional gas near uni- 
tarity [TU [34] . Equation ( 59 ) agrees with [34] and agrees 
with [T3] up to the same factor of 2/3 mentioned in the 
previous section. We obtain the high-frequency tail of 
the 3D bulk viscosity (60 1 for an arbitrary value of the 
scattering length. It vanishes at unitarity, in agreement 
with [14] . 



D. Viscosity sum rules 

By integrating over the frequency of the spectral func- 
tions, the asymptotic behavior of the correlation func- 
tions can be used to derive sum rules. On the one hand, 
such sum rules provide constraints that serve as a consis- 
tency check for approximate methods [34] . On the other 
hand, they are often more closely related to experimen- 
tal results since it can be difficult to probe the strictly 
asymptotic region of the response functions. The sim- 
plest way to obtain sum rules is to calculate the lowest 
momenta of the spectral functions. Additional sum rules 
can be deduced by introducing some form of UV cut- 
off that shifts the integration region to lower frequencies 

For completeness, we calculate the lowest moment of 
the bulk and the shear viscosity and rederive the results 
of Rcf. [T3] using the OPE. To obtain a finite expression 
for the shear viscosity sum rule, we subtract the high- 
frequency tail in Eq. ( 59 ) from the integrand. This yields 



do; 



C 



"■3D 
67TTO 



C (63) 



and 



1 1 9 97rm 



3D 



C 



lim lim 



m z uj 



ImG2(q, W ). (64) 



The shear viscosity sum rule ( 63 ) agrees with the result 
in [33] and agrees with [lj] up to the pref actors. 

The OPE alone is not sufficient to constrain the bulk 
viscosity sum rule ( [64] ) . As it was pointed out in Ref. [14] , 
this is due to an ordering ambiguity when taking the zero- 
momentum limit q — > and performing the u> integration 
(we take the limit after performing the integration) . The 
term in the second line of Eq. (64) is evaluated in the 



zero-momentum and zero-frequency limit where the OPE 
cannot provide an insight. It was evaluated explicitly 
in [14] using hydrodynamic expressions for the current 
response. Combining the result in [14) with Eq. (64) 
yields the sum rule 



dui £(w) 



1 



<9(a 3 D 



C) 



367rm 



da 3 u 



where the subscript s indicates that we take the deriva- 
tive of C with respect to at fixed entropy. 



IV. SUMMARY 

In this paper, we studied the current response func- 
tion of a Fermi gas and calculated the first terms of 
its operator-product expansion. We considered the limit 
uj, q — > oo away from the one-particle peak ui 7^ q 2 /2m 
with q 2 /mui and a _1 /V™ held fixed. We found that, to 
leading order, the response is proportional to the contact 
whose Wilson coefficient we determined. We calculated 
various response functions for a wide range of parameters: 
we derived the longitudinal and transverse components 
of the spin-symmetric and antisymmetric response func- 
tions in 3D and 2D for arbitrary values of the scattering 
length. 

Using those results, we were able to obtain the asymp- 
totic form of the dynamic structure factor in 3D and 2D 
and thus to generalize and to extend many of the pre- 
vious results for 3D systems at large scattering length 
[T2l [T4 , 22 . We rederived the high-frequency tails of the 
bulk and shear viscosity considered in [14J using the OPE 
and calculated their 2D counterparts. Furthermore, we 
could use the OPE to calculate some sum rules that con- 
strain the spectral viscosities in 3D. 

As illustrated by this work, the operator-product ex- 
pansion provides a powerful tool to derive universal re- 
lations. It offers a systematic way to analyze the short- 
time and distance structure of a theory, which can be 
fully characterized by performing calculations involving 
few-body states. 

Note added. — Recently, Ref. [35] appeared. The au- 
thors calculate the viscosity tails in 3D and obtain a re- 
sult in agreement with Eq. (59) and (60) of this paper. 
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Appendix A: MATRIX ELEMENTS FOR THE 
TWO-PARTICLE STATE 

This appendix gives some details of the calculation of 
the Feynman diagrams in Fig. [4] We begin by sketching 
the most important steps in the evaluation of the ma- 
trix elements between the two-particle state. Section [A~T] 
gives an example of such a calculation and shows how to 

(3) 

derive the contribution Cy„ to the Wilson coefficient of 

t j 

the contact operator [Eq. (37)1. The results for the re- 



(65) 

maining contributions Cff and are listed in Sec. 



A2 



We evaluate each diagrams using the Feynman rules 
summarized in Sec. |IIB] For all diagrams involving an in- 
ternal loop with loop momentum k and energy /cqj the fco 
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integration can be carried out trivially using the theorem 
of residues, which gives the expressions in Eqs. (35 1, (36 ), 



and (37). In the next step, we perform a partial fraction 



decomposition of the integrand with respect to k. Where 
necessary, we introduce Feynman parameters and shift 
the integration variable to eliminate the angular depen- 
dence in the denominator. Since the resulting expres- 
sion is rotationally invariant, we only need to retain the 
terms that are symmetric under k — > — k. Rotational 
invariance also dictates the tensorial structure of the in- 
tegrals. That way, we can express the integral as a sum 
containing (scalar) integrals of the type 



k" 



d d k 



(2n) d (k 2 - a 2 - is 



(Al) 



where a 2 is some constant (that typically depends on u>, 
q, and on a Feynman parameter). The value of those 
integrals can be determined by performing the integral 
for arbitrary a and 6=1 and differentiating the results 
with respect to a 2 . Of course, one can also analytically 
continue a — > ia and evaluate the resulting integral. It is 
convenient to introduce IR and UV cutoffs ij and A, re- 
spectively, when performing the calculations. It remains 
to calculate the integrals over the Feynman parameter. 
They take the form 



d.r 



{1-xf 



o ' [x(mu - q 2 /2) + x 2 q 2 /A] a ' 



(A2) 



If we define the variable z = q 2 /4muj and perform the 
transformations z — > z' = z/(l — z) and then x — > 1 — x, 
we can recast Eq. ( A2 ) in the form 



1 , x b (l-x)- a 
dx 



(mcu - q 2 /4) a J [1 - xz'] a 

1 r(i + 6)r(i-a) 

~ (mu-q 2 /4) a T(2-a + 6) 
x 2 F 1 (a,b + 1;2 - a + b; z') , 



(A3) 



where the left-hand side defines the Gaussian hypergeo- 
metric function 2-F1 |36j . The functions can be further 
simplified using standard identities such as |37) . We can 
rewrite the resulting expression in terms of elementary 
functions using any computer algebra system. In 2D, ad- 
ditional integrals may arise where the integrand depends 
logarithmically on x. They can be solved with the aid of 
a computer algebra system as well. 

While the evaluation of each Feynman diagram pro- 
ceeds as outlined above, the calculations can be very cum- 
bersome and will not be reproduced here. As an example, 
we consider only the simplest diagrams and calculate the 
matrix elements of the three diagrams in Fig. 4(c) in the 



next section. The results of the other diagrams are given 
in Sec. [O 



(3) 

1. The contribution C"y to the contact's Wilson 
coefficient 



To illustrate the procedure outlined in the previous 



section, we calculate the sum of the diagrams in Fig. 4(c 



If we apply the Feynman rules in Sec. II B to evaluate 



Fig. |4(c) we obtain 
(3) _ iA((J,q) 



C 



A- 



u-E, 



-1i 



dfc 
2tt 



d d k 

(27rp 



fc - E k 

i 



E, 



is k +uj - E k+ q 
-{2ki+q t 



IE 



IE 



3] + (w 



(A4) 



where we divided by the matrix clement of the contact 

(3) 

operator in order to obtain the contribution to the 



contact's Wilson coefficient. Equation (A4) can be fur 



ther evaluated using the theorem of residues 



r<(3) _ 



iA(u,q) 



u — E, 



-A~ 



1 



1 



x 



2E k lo - E'k+q - E k + is 



d d k 

(27r) d 

(2fc< + qi 



(A5) 



which coincides with the expression in Eq. (37). 



In the next step, we introduce Feynman parameters 
and shift the integration variable k — > k — xq/2 to obtain 
a rotationally invariant integrand. The term ~ fcj is odd 
under k — > — k and can be dropped. This yields 



C\f = -iA{uj, q) 



1 1 
mu> — q 2 /2 mA 
1 



dx(l — x) 



(2n) d {k 2 - (x(muj - q 2 /2) + x 2 q 2 /4) - is} 2 _ 
+ (u}->-w). (A6) 

In order to perform the k integration, we note the fol- 
lowing identity: 



d d k 1 

(27r) d k 2 - a 2 



IE 



Si(A+fa) d = 3 
^(iln^ + f) d = 2 



(A7) 

Differentiating this with respect to a 2 gives the desired 
result 



d d k 



1 



(2Tr) d [k 2 -a 2 - is} 2 



where we have to put a 2 = x(moj — q 2 /2) + x 2 q 2 /A. We 
consider the 3D case first. Inserting Eq. ( A8 ) in ( |A6[ ) and 
using (A3) gives 



Cf = -W(w,q) 



1 1 

muj — q 2 /2 mA 



13 



67T ^Jmuj - q 2 /A 



(A9) 



The hypergeometric function can be written in terms of 
elementary functions as follows: 



2-Fi ( — i 2: — : z' 
~ 1 2 2 



, 3mw / 277 I i mw + <jV ma; ~ <7 2 / 4 \ 

+ —3- V» - T / 4 In Y 27^ ■ 

<7 J V tow — g z /2 / 



(A10) 



Substituting this in Eq. (A9) yields the final result for 
the Wilson coefficient 



C 



(3) _ iA(u, q) 



4tt 2 



X 3D 



2(mu}-q 2 /2) q 2 
+ qy/ rnuj — q 2 /4" 1 2 



+ -^\Jmuj - q 2 /A 



muj muj 

— o — In 

q 3 mu! — q 2 /2 



q l q J 



(LU — > —Ld). 

(All) 



Similarly, we obtain for the 2D case: 



a 



(3) 



A(u, q) 



13 4tt 2 {m,uj-q 2 /2y 



In a 2 D v— (mw - <7 2 /4) 



1 2tow mu — q /A 



2/Al 2 



In 



tow — (? 2 /2 



(A12) 



2. List of contributions to the contact's Wilson 
coefficient 



For completeness, we state the contributions of the 
other diagrams in Fig. [4] to the Wilson coefficients with- 
out giving an explicit derivation. Figure 4(a) contributes 
a term 



C 



(i) 



1 



Airmq 2 \ (tow - q 2 /2) 
1 



^ -f/* (g 4 - 4tow, 2 + 6to 2 w 2 ) *^ In ^ + *V^W±\ 



tow — q 2 /2 



2-Kmq 2 



I 5— mu tow + qy/mu)- q 2 /A 

y/mu) - H In — 

q tow — q 1 2 



Sij ~ c«(q, w) iSD + cg>(q, w)|» + ( W -+ -w) 

47TTO J 27TTO 

(A13) 



in 3D. & is defined in Eq. |29| and cf) in Eq. In 2D, we find 



2™, , 7^4 



1 / 42g 2 TOw - 7q 
47TTO mw - q 2 j2 \ 48(mw - g 2 /2) 2 



tow) 2 4tow mw-q 2 /4 — — 1 . 

In Vtt: - 2 In a 2D V - ( mw - 9 / 4 ) ) 



g 2 mw — q 2 /2 



2 ^2 / Yytiuj Wq 2 2ttl 2 uj 2 

7 7 27^71 "1 77 2 h ( 3muj ~ ? ) ln ° 2D V-( muj - <7 2 /4) 

47TTO (tow — q I2y \ 4 24 g- 2 v ' 



tow(8to 2 w 2 — 12towq 2 + 3q 4 ) tow — q 2 /4\ 



24 

q 4 muj — q 2 /2 J q 2 

Similarly, we obtain for the sum of the diagrams in Fig. |4(b)| 

•2 



+ (w -> -w). 



(2) 



1 I 1 ^274 _|_ TOW / 

47TTO |_ q 2 ' q 3 \ 

3(tow) 2 q 2 



(tow) 2 



TOW 



hi 



tow + q^Jmui — Q 2 /4 



tow — g 2 /2 



47TTO 2 Wq 
„■„-! 



Ill 



tow + q^/mw — (? 2 /4 3tow 



tow — <? 2 /2 



\/ tow — 9 2 /4 



gigj 

/t2 



3D 



-«f + (w^-w) 



in 3D, and 



C (2) = 

lJ 47TTO 2 W TOW — ? 2 /2 



In 



87TTO 2 W TOW — <7 2 /2 

to 2 w 2 — Amujq 2 + q 4 tow — g 2 /4 
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In — - In a 2D \/-{mu - q 2 /4) 

tow — q z / 2 q z 



A 



1 



27TTO 2 w V y/-(muj - <? 2 /4) 



hi 



tow — q 2 /4 
muj — q 2 /2 



Sij + (w — > — w) 



(A14) 



(A15) 
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(A16) 
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in 2D. 
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